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We use extensive Monte Carlo transfer matrix calculations on infinite 
strips of widths L up to 30 lattice spacing and a finite-size scaling analysis 
to obtain critical exponents and conformal anomaly number c for the two- 
dimensional XY-Ising model. This model is expected to describe the critical 
behavior of a class of systems with simultaneous U(l) and Z2 symmetries of 
which the fully frustrated XY model is a special case. The effective values 
obtained for c show a significant decrease with L at different points along the 
line where the transition to the ordered phase takes place in a single transition. 
Extrapolations based on power-law corrections give values consistent with 
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c = 3/2 although larger values can not be ruled out. Critical exponents 
are obtained more accurately and are consistent with previous Monte Carlo 
simulations suggesting new critical behavior and with recent calculations for 
the frustrated XY model. 
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I. INTRODUCTION 



Recently, the critical behavior of the two-dimensional XV-Ising model, consisting of 
XY and Ising models coupled through their energy densities, has been studied in some 
detailEH. The model is expected to describe the critical behavior of a class of systems with 
U(l) and Z 2 symmetries which includes, for example, two-dimensional fully frustrated XY 
(FFXY) modelsii'iii'ii'0'0iii'lli'0, or alternatively, two-dimensional arrays of Josephson 
junctions^!, one-dimensional ladders of Josephson junctions with charging effects^, helical 
XY mo delH and some surface solid-on-solid modelsE3H3. The Xl^-Ising model is also of 
great theoretical interest in its own right as the phase diagram obtained by Monte Carlo 
simulations gave rise to interesting and unusual critical behaviorS. Recent work by Knops et 
alM has further justified the relation between XF-Ising and FFXF models by showing that 
the phase coupling across chiral domains in the FFXY model is irrelevant at criticality. In 
the subspace of parameters of the model where the XY and Ising coupling constants have 
the same magnitude, separate XY and Ising transitions, first-order transitions and a critical 
line with simultaneous XY and Ising ordering were found. 

The numerical study revealed that starting at the branch point where separate XY and 
Ising transitions merge, the line of single transitions has a segment of continuous transitions 
which eventually become first order as one moves away from the branch point. Along the 
segment of continuous transitions, the critical exponents associated with the Ising-like order 
parameter were found to be significantly different from the pure Ising values and, in fact, 
appeared to be non-universal, varying along the line. Besides by critical exponents, this 
critical line was also characterized by its central charge, or conformal anomaly number c. 
The central charge was estimated from the finite-size scaling of the free energy of infinite 
strips at criticality, obtained from Monte Carlo Transfer Matrix calculations. The results 
obtained from strips of width up to L = 12 lattice spacing were rather surprising: the central 
charge appears to increase continuously along this line, from c ~ 1.5 close to the branch 
point to c ~ 2 near the tricritical point. 
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Similar calculations00'El of the critical exponents and central charge for the FFXY 
model were consistent with these results. Although models with varying c are well known, 
as for example the g-state Potts and 0(n) models with a continuously varying number of 
states q and n, the behavior for XF-Ising model is rather unexpected since, contrary to 
the previous models, the transfer matrix can be chosen symmetric and along the critical 
line a parameter is changing that does not affect the symmetry. The question then arises 
if this behavior is a real effect or an artifact due to limited strip widths. In view of the 
relation between the XF-Ising and FFXY models, the answer to this question may also 
give some insight into the behavior of the central charge for FFXY models@0lil0. Also, 
it is important to have improved estimates for the critical exponents in order to be more 
certain about the non-Ising nature of the critical behavior along the line of single transitions. 

In this work we report the results of extensive Monte Carlo transfer matrix calculations, 
using infinite strips of widths L up to 30 lattice spacing, aimed to resolve some of the 
issues raised by previous calculations on the AF-Ising model. Rather than attempting to 
evaluate critical exponents and central charge at several different points along the line of 
single transitions to check if these quantities do change or remain constant along the critical 
line, we have concentrated on a couple of points but performed extensive calculations for 
large L and used variance reduction techniques to decrease the statistical errors. The results 
for the effective value of c show a significant decrease with increasing L, indicating that they 
even the extrapolated estimates have not yet reached their asymptotic values for L = 30, 
the largest strip width considered. Extrapolation suggests values not inconsistent with 
c = 3/2. However, on purely numerical grounds, we can not rule out the possibility of a 
larger value or even a varying c along the line. Our results for the central charge suggest 
that the recent estimates of this quantity, c ~ 1.6 for the related FFXY mo delsfflffl 
are likely to be subject to similar systematic errors due to slowly decaying corrections to 
scaling and the asymptotic value is in fact consistent with c = 3/2. The critical exponents 
associated with Ising-like order parameter are obtained more accurately, although there are 
some puzling inconsistencies. The exponents are found to be significantly different from the 



pure Ising values but consistent with the previous Monte Carlo simulations which suggested 
new critical behavioii and recent estimates for the FFXY model using; Monte Carlcll and 
exact numerical transfer matrix calculations^. 

The paper is organized as follows. In Sec. II, we define the model and briefly review 
the main features of its phase diagram, indicating the locations near the phase boundary 
where the Monte Carlo transfer matrix calculations were performed. In Sec. Ill we provide 
details on the Monte Carlo transfer matrix method and the implementation of the variance 
reduction techniques. In Sec. IV, a finite-size scaling analysis of the interfacial free energy 
is used to extract critical quantities. In Sec. V, we present the numerical results for critical 
couplings, exponents and central charge and in Sec. VI we discuss and compare these results 
with previous calculations. Finally, Sec. VI is devoted to the conclusions and final remarks. 

II. MODEL AND PHASE DIAGRAM 

The W-Ising model is defined by the following HamiltonianB! 

= - Yl K A + Bs i- s Mi ■ n i + Cs iSj ], (1) 

K1 <ij> 

where s — ±1 is an Ising spin and n = (cos 8, sin 9) is a two-component unit vector, is an 
XY spin. The model can be regarded as the infinite coupling limit, h — > oo, of two XY 
mo delsiii coupled by a term of the form hcos2(6 1 — 6 2 ) and has a rich phase diagram in 
the A, B plane that depends strongly on the value of C . The model with A ^ B is relevant 
for the anisotropic frustrated XY model! and anti-ferromagnetic restricted solid-on-solid 
mo del0. 

In this work we will be concerned with the critical behavior of the W-Ising model of 
Eq. ([]]), defined on a square lattice, in the subspace A = B, 

T7f = -J2 W + S i S i) n i ■ n J + Cs i s i], ( 2 ) 

which is relevant for the isotropic frustrated XY model or its one- dimensional quantum 
version@. The phase diagram obtained by Monte Carlo simulations is shown in Fig. [I] and 



consists of three branches joining at P, in the ferromagnetic region A > 0, A + C > 0. 
One of the branches corresponds to a single transition with simultaneous loss of XY and 
Ising order and the other two to separate Kosterlitz-Thouless (KT) and Ising transitions. 
An important feature of the phase diagram is that there is no phase with Ising disorder and 
XY order thus indicating that Ising disorder induces also XY disorder in this model. This 
is related to the special symmetry under the transformation 

uj -> sjiij (3) 

which holds if A = B, since XY spins are not coupled across an Ising domain wall where 
SiSj + 1=0. The behavior of FFXY model coresponds to the behavior of this model along a 
particular path through this phase diagram. The available numerical results for the standard 
FFXY modeifflS are consistent with a single transition but generalized versions could 
correspond to a path through the double transition region. In fact, a Coulomb-gas represen- 
tation of the FFXY model with an additional coupling between nearest neighbor vortices 
has a phase diagram with identical structured. In the one-dimensional quantum version of 
the frustrated XY modefl, related to the zero-temperature transition of Josephson-j unction 
ladders, double or single transitions will result depending on the ratio between inter-chain 
and intra-chain couplings. In the Monte Carlo simulationsi, the critical line PT in Fig. [I] 
appears to be non-universal as the critical exponents associated with the Ising order param- 
eter were found to vary systematically along this line. In addition, a preliminary evaluation 
of the central charged c using data for the free energy of infinite strips obtained from Monte 
Carlo transfer matrix appeared to indicate that c varies from c ~ 1.5 near P to c = 2 near 
T. These results for the central charge were based on strips of width L up to 12 lattice 
spacings. However, this range of L and the numerical noise in the data does not allow one 
to extrapolate to the large L limit and these estimates are thus subject to systematic errors. 
To obtain accurate estimates it is necessary to perform calculations for larger L and also to 
reduce the errors. Rather than attempting to evaluate c at several different points along the 
line PT in order to check if this quantity changes or remains constant along the line, we have 
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concentrated our attention at a few points but performed extensive calculations for large L 
and used variance reduction techniques to decrease the errors. The calculations discussed 
in the following sections were performed primarily along the cuts through the critical line 
as indicated in Fig. [[]. 

III. MONTE CARLO TRANSFER MATRIX 

Estimates of the free energy density per lattice site was computed using the Monte Carlo 
transfer matrix method. We give a brief summary of the essentials of this method and refer 
the reader to Ref. ^l] for more details. 

Helical boundary conditions are convenient for these computations. In this case the 
transfer matrix can be chosen to be a sparse matrix for the case where one matrix multipli- 
cation corresponds to the addition of one surface site to the lattice, as illustrated in Fig. ||. 
The sparseness follows from the fact that from any given configuration of surface sites only 
those new configurations can be reached that differ at that newly added site only. We used 
a transfer matrix defined by 

T(s u . . . , s L ; mi, . . . , m L \ t x , . . . , t L ; m, . . . , n L ) = (4) 

e (An L -i-n L +Bn L - 1 -n L t L - 1 t L +Ct L - 1 t L +An L -mi+Bn L -m L s L s 1 +Cs L - 1 s L ) j^L ^ (j(m k — n k ^S s j 

The statistical variance of transfer matrix Monte Carlo computations is proportional to 
the variance of the quantity 

^ = ^2 d/bfrnj T(si, ...,s L ; /bfm 1; . . . , m L 1 1 1} . . . , t L ; m, . . . , n L ). (5) 

si J 

The variance can be reduced by applying the transfer matrix algorithm to a similarity 
transform T of the transfer matrix T defined in Eq. (|5]). The transformation requires an 
optimized trial eigenvector ip^ and is defined as follows: 

f({s}; {m} | {t}; {n}) = ^ T ({,}; {m})T({,}; {m} | {t}; {n})/^ T ({t}; {n}) (6) 

where {x} = x±, . . . , xl- 
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In the ideal case, when ip? is an exact left eigenvector, the local eigenvalue /t, defined 
by Eq. (|^) with T replaced by T, is a constant - an eigenvalue of the transfer matrix. In 
practice, the better the quality of the trial function, the smaller the statistical noise in the 
Monte Carlo estimates of the transfer matrix eigenvalues. 

For the design of good trial states it is helpful to realize that the dominant left eigenvector 
is proportional to the conditional partition function of a semi-infinite lattice, extending to 
infinity towards the left, as indicated in Fig. ^ conditional on the microscopic state of the 
surface. 

Our computations used the following form for the trial vectors: 

V>t(si, • • • , s L ; m, . . . , n L ) = exp ^ (A itj n; • + B itj n; • + C^s^) j . (7) 

Here the parameters A+j, By, and CV,- are variational parameters, which are chosen so as 
to minimize the variance of ft, as described in detail in Ref. [21]. In the expression (0) the 
asterisk indicates that the sum over the pairs of surface sites labeled i and j is truncated, 
as is required for for computational efficiency To truncate in a way that allows systematic 
improvement of the quality of the trial function, it is necessary to guess for which pairs of 
sites % and j the interaction parameters A it j in Eq. (|7]) have the largest magnitudes, and 
similarly for B{j, and Cy . Obviously, interaction strengths will decay with distance, but 
owing to the helical boundary conditions and the surface defect, the geometrical distance is 
not quite correct. Instead, a distance can be defined between surface sites i and j of the 
semi-infinite strip (illustrated in Fig. as the length of the shortest path that: (a) connects 
sites i and j; (b) passes only through bulk sites (indicated by full circles in Fig. H); and (c) 
travels along the edges of the square lattice. 

The reason for excluding surface sites from the path is that the correlations described 
by the interaction parameters Ay, etc. are mediated only via bulk sites, since those are the 
only ones that contain variables that are not frozen in the conditional partition function. 
Fig. ^| shows a path of length three. Owing to the presence of the surface defect no two 
paths are strictly equivalent and, since the surface interactions can be regarded as functions 
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of the minimal path defined above, all parameters have to be assumed different. However, 
the transformed transfer matrix T depends only on the ratio of the values of two trial 
states shifted by a single lattice unit along the surface. By artificially imposing translational 
invariance on the interaction parameters, one can produce cancellations in the computation. 
This reduces the number of arithmetic operations from order L (in the absence of translation 
symmetry) to a number of the order of the maximum path length at which the interaction 
are truncated. 

Suppose that interactions in the trial function are truncated at path length /, measured 
in units of the lattice spacing, then the following compromise seems to work satisfactorily: 
give those interaction parameters the same values that are (a) farther away from the defect 
than Z, and (b) would be equivalent by translation symmetry in the case of simple periodic 
boundary conditions. In particular, this means that all interactions associated with paths 
that cross the defect are allowed to be different in the computation. It should be noted that 
this approximation can only be improved to a point: as soon as many-body interactions 
appear that are of greater strength than pair interactions included in the trial function, 
ignoring the many-body interactions makes it impossible and pointless to determine the 
two-body interactions. 

We are interested only in studying the behavior of systems with B = C, but the twisted 
boundary conditions force us also to consider the case where B = —C, which is obtained by 
inverting either the Ising or the XY variables on one sub-lattice. In all of these cases we used 
trial functions in which the corresponding relation was maintained between the interaction 
parameters appearing in the trial vector, Eq. (|7|), i.e., B^ = Cy if i and j belong to the 
same sub-lattice and Bij = —Cij otherwise or if twisted boundary conditions are used. 

As a final comment we mention that by using the variance reduction scheme mentioned 
above the Monte Carlo calculation can be accelerated roughly by a factor of two hundred^. 
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IV. FINITE-SIZE SCALING 



To locate the critical couplings and determine the critical exponents we will do a finite- 
size scaling analysis of interfacial free energies. Since the model contains both XY and Ising 
variables, there are in principle two types of interfacial free energies that can be determined 
by a suitable choice of the boundary conditions. If a twist in the Ising variables is imposed 
by anti-periodic boundary conditions, a domain wall is forced along the infinite strip and the 
associated interfacial free energy can be obtained from the difference per surface unit the 
free energies of systems with periodic and anti-periodic boundary conditions. On the other 
hand, if the same procedure is followed for the XY variables, a smooth phase twist of tt is 
forced across the infinite strip. The transfer matrix calculations are done for an L x oo strip 
with L even and helical boundary conditions. With this set-up, it is simple to introduce 
independent twists in the Ising and XY degrees of freedom. 

The interfacial free energy of an Ising domain wall of length L along the strip is given by 

AFi = L 2 [f(A, A, C) - f(A, -A, -C)], (8) 

where f(A, B, C) is the free energy per site of the Xy-Ising model of Eq. ([J). The parameters 
A and C are chosen so that the ground state is ferromagnetic, A > and A + C > 0, so 
that taking B — > —B = —A induces a domain wall between the two anti- ferromagnetic Ising 
ground states. Similarly, a twist of 7r in the XY degrees of freedom is induced by A — > —A 
and B — > —B so that 

AF XY = L 2 [f(A, A, C) - f(-A, -A, C)], (9) 
and the helicity modulus 7 is given by 

7 = 2AFxy/7T 2 . (10) 

At a conventional second-order transition, the interracial free energy has the simple 
scaling form 
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AF(A, C; L) = A(L yT t), (11) 

where A(u) is a scaling function and t(A, C) is the non-linear scaling field measuring the 
distance from the critical point; the thermal scaling exponent is related by v = l/yr to 
the exponent is, which describes the divergence of the correlation length at criticality. In our 
analysis, we fix one of of the parameters, A or C, and expand t to quadratic order, i.e., for 
fixed C we have t = A — A C (C) + k[A — A C (C)} 2 and similarly when A is kept fixed. In the 
vicinity of the critical coupling t — 0, a standard finite size scaling expansion in u = tL VT 
for u ~ yields 

AF(A, C; L) = a Q + a x u + a 2 u 2 H . (12) 

With our convention, u is positive in the ordered phase, AF will increase with L for u > 0, 
decrease for u < and be a constant at u = for sufficiently large L so that corrections to 
scaling have become negligible. 

Sufficiently close to u = 0, Eq. QT2| ) can be used to obtain accurate estimates of the 
critical exponent i/t (or equivalently is) and the critical values A and C. The expansion is 
truncated at some high order (u 5 in some cases). A critical dimension x^ d ' of a disorder 
operator can be obtained from the constant ao via 

«« = £. (13) 

Z7T 

The critical dimension a; describes the decay with distance r at criticality of the two-point 
correlation function g(r) of an operator determined by the choice of boundary conditions: 
g{r) oc r~ 2x(d) . The scaling exponent = 2 — describes the behavior of this operator 
under scaling. 

As mentioned above, we consider two kinds of anti-periodic boundary conditions. Sub- 
scripts will be used to distinguish the exponents and of the associated operators. In 
the case of anti-periodic boundary conditions in the XY variables, the conjugate operator 
is a vortex of strength ~, measured in units 2n. Such an operator is represented as the end 
point of a path on the dual lattice: XY bonds crossing this path have their interactions 
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changed from A to —A. Because of the symmetry of the model under the transformation 
given in Eq. (|3]) this operator is equivalent to one where the B is changed to —B. The 
exponents of the ^-vortex will be denoted by x\® i and y^L i . The operator corresponding 
to the case of antisymmetric boundary conditions for the Ising variables is the standard Ising 
disorder operator. The exponents for this case with be denoted and y^. For self-dual 
models or models for which a renormalization mapping to the Gaussian model exists, the 
disorder operators can be related to order operators, but we cannot derive either of those 
properties for the XV-Ismg model. 

The critical exponents were estimated by making a scaling plot of AF in which the 
parameters were estimated by a constrained least-square fits with the critical couplings 
fixed at their most reliable estimates, i.e., those obtained by extrapolation from the Ising 
domain wall data. Unfortunately, the discontinuity in the helicity modulus 7 is not accessible 
by similar finite size scaling considerations since the discontinuity in 7 is defined in the 
thermodynamic limit L —>■ 00 and 

A 7 = 2[.A(oo) - ^(-oo)]/tt 2 . (14) 

A rough estimate from Fig. |] for ^4(±oo) gives A7 « 1.3 which is about double the value 
2/tt of the two-dimensional XY model. This estimate is not very reliable but we can say 
with a considerable degree of confidence that A7 is much larger than 2/tt in this coupled 
XF-Ising model and in the FFXY model. 

In addition to critical exponents, another important quantity which provides information 
on the nature of the critical behavior is the central charge c of the conformal invariance. 
This quantity can be obtained from the amplitude of the singular part of the free energy 
per site@, at criticality, in the infinite strip by 

f(A c ,C c ,L)^f 00 + — (15) 

for sufficiently large L, where denotes the regular contribution to the free energy at 
the critical point. The central charge classifies the possible conformally invariant critical 
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behaviors. For example, for the pure Ising model, c = 1/2, and along the critical line 
of the XY model c = 1. Although c is only defined at criticality, Eq. (|T5| ) can be used 
to define a size and coupling dependent effective central charge c(A, C, L) away from the 
critical point. If this quantity is now identified as the function c(g) defined in the c-theorem 
of Zamolodchikovil, where g stands for a set of coupling constants, this quantity should 
have a well-defined behavior near criticality since, according to the c-theorem, c(g) is a 
monotonically decreasing function under a renormalization group transformation and reaches 
a constant value equal to the central charge at the fixed point. An interesting consequence 
of this identification is that c(A, C, L) should have a maximum near the critical line of single 
transitions in the AF-Ising model with a lower bound c(A, C, L) > 3/2 and away from 
the critical line should converge to either c(A, C,L) = 1 in the XY ordered phase or to 
c(A, C, L) = in the remaining phases. Our calculations are consistent with this behavior 
but we found that the maximum in c(A, C, L) does not provide an accurate procedure to 
locate the critical couplings since it is rather flat within a wide range of couplings near 
u = 0. To obtain an estimate of the central charge c at criticality we first accurately locate 
the critical couplings using the non-linear fitting of Eq. ( [12] ) and extract a size- dependent 
c(L) from the singular part of the free energy in Eq. (|T5|), which is subsequently extrapolated 
to L — > oo. 



V. ESTIMATES OF CRITICAL POINTS, EXPONENTS AND CONFORMAL 

ANOMALIES 

We computed eigenvalues of the transfer matrix for various points along the critical curve 
and used these to extract estimates for the central charge. In two cases we recomputed the 
critical points themselves from a scaling analysis of the interface free energy and helicity 
modulus. We start our discussion with the latter. Fig. |^ is a scaling plot of the Ising 
interface free energy as a function of A at fixed C = 0.2885. Fig. |] is the same for XY 
interface, obtained by choosing boundary conditions that induce a twist of ir in the XY 
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variables. Figs. || and ^| are analogous plots for the case A = 2 with varying C. The scaling 
plots for the systems with anti-periodic boundary conditions in the Ising variables do not 
display statistically significant deviations from the scaling hypothesis. However, this is not 
the case for the scaling plots for systems with a twist in the XY variables as shown most 
clearly by Fig. In fact, significant changes are ,e.g. in the critical amplitude, are observed 
in the "scaling plot" if smaller system sizes are omitted from the fit. 

In all cases, there are strong corrections to scaling for small systems. This is demonstrated 
in Figs. [7] and || plots of the estimated effective critical couplings versus inverse system size. 
The effective coupling associated with size L was obtained by a least-squares fit to system 
sizes including sizes L and up. 

By extrapolation assuming overly conservative 1/L 2 corrections, we obtain the following 
estimated critical points: A = 1.0014 (Ising twist) and A = 1.0025 (XY twist) at C — 
0.2885, where the first of these two is presumably the more reliable one. Similarly, for A = 2 
the results are C = 1.318 (Ising twist) and C = 1.316 (XY twist). Our estimates for the 
critical exponents are summarized in Table Q; the plots in Figs. [7| and |S] may serve to provide 
rough error bars. 

Finally, we estimated the conformal anomaly c along the critical line using Eq. (JTSJ) 
and taking 2, 3 or 4 consecutive system sizes. This defines an effective c(L) at L^ 1 = 
(L~J n +L~^ x )/2, where L min and I/ max are the smallest and largest size used in the fit. Critical 
points not mentioned above, viz. (A = 3, C = —2.3250) and (A = 0.6, C = 0.1520), were 
taken from the estimates provided in Refs. p] and 0. The results are summarized in Fig. [13. 



VI. DISCUSSION 

The results obtained from the Ising interface free energy, summarized in Table [I], seem 
sufficiently accurate to exclude pure Ising critical exponents (yx = 1 and yf^ = 15/8) for 
the point Awl. Our numerical results agree with those for the 19- vertex model obtained 
by Knops et ai.@, who find = 1.23(3) and y^ = 1.80(1). Within the sizable uncertainties 
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in the estimates of the same exponents for A = 2, we find no evidence for variation of these 
exponents along the critical line. The results obtained for the thermal exponent yx are 
consistent with those from direct Monte Carlo simulations§ of the AF-Ising model: 1.19(4) 
for A — 1 and yx = 1.18(4) for A = 2. We note, however, that the latter computations 
indicate variation along the critical line of the scaling exponent of the order parameter, an 
exponent which was not computed in the present transfer matrix Monte Carlo approach. 
Lee et al. found: y[ o) = 1.85(2) for A = 1 and y[ o) = 1.80(2). 

There is a serious internal inconsistency in our estimates for yx as obtained from the 
the Ising interface and those obtained from the XY interface. Although Figs. || and [LC 
display strong corrections to scaling, there is no indication that the two ways of computing 
this thermal exponent will become consistent in the limit of large systems. This calls in 
question the validity of the basic assumption of scaling theory, viz., that there is a single 
divergent length scale in this system as the critical point is approached along a temperature- 
like direction. 

The results for the critical exponents yx and yj for the XF-Ising model are consistent 
with similar Monte Carlo transfer matrix calculations for the FFXY model on a square 
lattice^ (y T = 1.25(6), yf ] = 1.81(2) ) and its one- dimensional quantum versionlll (y T = 
1.24(6), yf ] = 1.77(2)), although the strip widths are much larger here and the accuracy 
much better. Estimates of the critical exponent yx for the FFXY model obtained from 
Monte Carlo simulations^, 1.21(3) (square lattice) and 1.18(3) (triangular lattice) are also in 
good agreement with the result for the XY-Ising model and seem to support an XF-Ising 
universality class for these systems. 

The exponent yf 1 = 1.85 for the Xl^-Ising model obtained from Monte Carlo 
simulations^ is significantly larger than the result for y^f 1 in Table I. This discrepancy is 
also observed in the results for the FFXY modeJ^B and can in part be attributed to cor- 
rections due to the effective free boundary conditions for the XY degrees of freedom at 
criticality. As argued in the context of the FFXY mo since the XY degrees of free- 
dom are uncoupled across an Ising domain wall, the XY variables should be regarded as 
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having free boundary conditions instead of periodic ones. This results in a correction to 
the estimate of — > — 1/16 which seems to improve the agreement between y\°^ and 
y^ although, as mentioned in Sec. IV, the precise relation between these exponents is not 
known. 

The results for the exponents in Table I and the scaling plots for A(u) for the interface 
free energies are based on the naive assumption that the XY and Ising correlation lengths 
behave as £ a ~ t~ l ^ Va with y a determined by independent best fits for the Ising and XY 
interfacial free energies. Such a procedure would certainly be incorrect if the XY and Ising 
variables were decoupled as then log£xy ~ t~ l l 2 . In the present case, these degrees are 
strongly coupled and the appropriate scaling form is not known and the XY degrees of 
freedom are probably subject to large slowly decaying corrections-to-scaling making the 
analysis of the data fraught with difficulty and uncertainty. A detailed analysis of the data 
for the XY twist free energy shows that the estimates of the finite-size scaling parameters Oj 
of Eq.(12) are not stable and depend on the number of data points included in the fit. For 
this reason, the scaling plots of A{u) are somewhat misleading and a naive use of Eq.(13) to 
estimate Xxy yields the L = 4 value for = 2 — x^ of Figs. (11, 12). The reason behind 
this is that the small L data has the lowest x 2 an d is most heavily weighted in the scaling 
plots of Figs. 4 and 6. It is amusing to note that the use of these estimates together with the 
with the relation, valid for a self dual Gaussian model 14 , ^xr,i/2 a; xy,i = 1/16, gives results 
in agreement with those of Knops et al 14 . However, we consider this to be fortuitous and 
not to be taken seriously. Another difficulty with analyzing numerical data for the XY twist 
free energy is that there must be a crossover to a low temperature Gaussian fixed line when 
u » as the low temperature phase must be just a XY model when there is long-range 
Ising order. 

The well-known difficulties of analyzing numerical data for the helicity modulus in this 
system are compounded by this cross-over so it is not surprizing that we are unable to make 
definitive statements about the critical exponents for the XY variables. One might try using 
a dual roughening representation, but there are negative Boltzmann weights at the critical 
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point in the dual representation which will lead to some difficulties. Despite being able to go 
to relatively large strip widths of L = 30 we are unable to reach definite conclusions about 
the critical behavior of the XY degrees of freedom except to say that our simple ansatz 
for the scaling of the XY twist is inadequate and corrections to scaling should be included 
in the analysis, but we do not know the form these should take. Also, we are unable to 
estimate the discontinuity wA'y in the helicity modulus except to say that, at the critical 
point 7T7 m 1.1 and, at T~, irA'y « 4 which we believe to be a fairly realistic estimate. 

We now consider the results for the central charge in Fig. The results for A = 
0.6 correspond to the branch point in Fig. [I] as estimated from Monte Carlo simulations!. 
Convergence is very poor in this case. The effective c starts at c = 1.5 for small systems, 
decreases very slowly for intermediate systems and them decreases rapidly for the largest 
system sizes. It is not possible to estimate the asymptotic value for this case. In fact, 
this behavior suggests inaccuracy in the estimate of the critical point. The other curves in 
Fig. [1] correspond to different points along the line of single transitions. Again, corrections to 
scaling are decaying too slowly as a function of system size to allow an accurate estimate of c 
in the large L limit. However, c = 3/2 along the line is not inconsistent with the data. This 
is shown in Table |TJ where the central charge is estimated assuming power-law corrections 
of the form a/L 3 ~ s + /3/L 4 ~ s the leading correction to the free energy, ttc/6L 2 , in Eq. [15|. 
We chose s = 0.2 so that we could simultaneously fit the results for A m 1 and A = 2. 
This value, c = 3/2, would be the expected one if the critical behavior along the single line 
could be described as a superposition of critical Ising and XY mo dels§. Our results for 
the critical exponents y t and however are not consistent with this hypothesis and suggest 
that the coupling between the Ising an XY degrees of freedom is vital. The results of the 
extrapolations should also be viewed with caution since they are not completely justified. 
There could be other corrections as exp(— aL) or InL/L but due to the noise in the data 
any attempt to include such terms in the extrapolation is pointless. 
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VII. CONCLUSIONS 



We have obtained critical exponents and the central charge for the XF-Ising model using 
Monte Carlo transfer matrix calculations on infinite strips of widths L as large as 30 lattice 
spacings. The results for c show a significant decrease with increasing L but converge very 
slowly to an asymptotic value. An extrapolation procedure indicates that these values are 
not inconsistent with c = 3/2. However, on purely numerical grounds, we can not rule 
out the possibility of a larger value or even c varying along the line of single transitions. 
Our results for the central charge suggest that the recent estimates of this quantity for the 
related FFXY models are likely to be subject to similar systematic errors due to slowly 
decaying corrections-to-scaling and the asymptotic value is consistent with c = 3/2. The 
critical exponents associated with Ising-like order parameter are obtained more accurately 
and are found to be significantly different from the pure Ising values but are consistent with 
previous Monte Carlo simulations which suggested new critical behavior and also with recent 
estimates for the FFXY model using Monte Carlo and exact transfer matrix calculations. 
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Table Captions 



1. Critical exponents associated with the variables for which the boundary conditions 
were twisted. 

2. Estimates of the central charge c assuming the free energy per site to be of the form 
foo + + 77^ + 1^- Fits were made using data for L, L + 2, . . . , 30. For L > 10, (3 
was fixed at the value obtained from the L = 4 fit. In all cases the normalized x 2 was 
of order unity. 

Figure Captions 

1. Phase diagram of the XF-Ising modeffl. Solid and dotted lines indicate continuous 
and first-order transitions, respectively. Filled circles indicate the locations where the 
present calculations were performed. 

2. Left: graphical representation of the conditional partition function of a semi-infinite 
strip with helical boundary conditions, i.e., the left eigenvector of the transfer matrix, 
which is shown on the right. In the lattice on the left, open circles represent sites with 
variables tj and (i — 1, . . . , L) that specify the surface configuration upon which the 
conditional partition function depends. The full circles represent sites with variables 
that have been summed over. Right: graphical representation of the transfer matrix. 
The variables associated with the circles make up the left index of the matrix; the dots 
go with the right index. Coincidence of a circle and a dot produces a product of two 
6-f unctions. 

3. Scaling plot of the interfacial free energy with Ising-twisted boundary conditions; A w 1 
is varied at constant C = 0.2885. 

4. Scaling plot of the interfacial free energy with Xl^-twisted boundary conditions; A pa 1 
is varied at constant C = 0.2885. 
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5. Scaling plot of the interfacial free energy with Ising-twisted boundary conditions; C ~ 
1.32 is varied at constant A = 2. 

6. Scaling plot of the interfacial free energy with XY- twisted boundary conditions; C 
1.32 is varied at constant A = 2. 

7. Effective critical couplings A c vs. 1/L and the results of extrapolation to L = oo at 
C = 0.2885 for both boundary conditions. 

8. Effective critical coupling C c vs. 1/L and the results of extrapolation to L = oo at 
A = 2 for both boundary conditions. 

9. Effective yr vs. 1/L for critical point at A w 1. 

10. Effective vs. 1/L for critical point at A = 2. 

11. Effective vs. 1/L for critical point at A w 1. 

12. Effective vs. 1/L for critical point at A = 2. 

13. Effective conformal anomaly vs inverse system size 1/L for various values of A. 
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TABLES 

TABLE I. Critical exponents associated with the variables for which the boundary conditions 
were twisted. 
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TABLE II. Estimates of the central charge c assuming the free energy per site to be of the 
form /oo + # + X§-f + tJtf- Fits were made using data for L, L + 2, . . . , 30. For L > 10, /3 was 
fixed at the value obtained from the L = 4 fit. In all cases the normalized x 2 was of order unity. 
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FIGURES 
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FIG. 1. Phase diagram of the XY-Ising modelfeM Solid and dotted lines indicate continuous 
and first-order transitions, respectively. Filled circles indicate the locations where the present 
calculations were performed. 
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FIG. 2. Left: graphical representation of the conditional partition function of a semi- infinite 
strip with helical boundary conditions, i.e., the left eigenvector of the transfer matrix, which is 
shown on the right. In the lattice on the left, open circles represent sites with variables t{ and nj 
(i = 1, . . . , L) that specify the surface configuration upon which the conditional partition function 
depends. The full circles represent sites with variables that have been summed over. Right: 
graphical representation of the transfer matrix. The variables associated with the circles make up 
the left index of the matrix; the dots go with the right index. Coincidence of a circle and a dot 
produces a product of two 5-functions. 
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FIG. 3. Scaling plot of the interfacial free energy with Ising-twisted boundary conditions; A « 1 
is varied at constant C = 0.2885. 
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FIG. 4. Scaling plot of the interfacial free energy with X^-twisted boundary conditions; i^l 
is varied at constant C = 0.2885. 
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FIG. 5. Scaling plot of the interfacial free energy with Ising-twisted boundary conditions; 
C ?3 1.32 is varied at constant A = 2. 
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FIG. 6. Scaling plot of the interfacial free energy with Xl^-twisted boundary conditions; 
C ps 1.32 is varied at constant A = 2. 
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FIG. 7. Effective critical couplings A c vs. 1/L and the results of extrapolation to L = oo at 
C = 0.2885 for both boundary conditions. 
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FIG. 8. Effective critical coupling C c vs. 1/L and the results of extrapolation to L = oo at 
A = 2 for both boundary conditions. 
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FIG. 9. Effective t/T vs. 1/L for critical point at A « 1. 
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FIG. 10. Effective yx vs. 1/L for critical point at A = 2. 
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FIG. 11. Effective vs. 1/L for critical point at A « 1. 



30 



1.9 
1.85 

1.8' 
1.75 

1.7 
1.65 

1.6 
1.55 

1.5 





1 1 

<9> <$><$> ^ 




i i 


i 


- 






Ising twist kOm 
XY twist :+■ ■: 


- 




+ ++ + 


+ 


+ 


+ - 




i i 






i 



0.05 0.1 0.15 

L- 1 



0.2 



0.25 



FIG. 12. Effective y (d) vs. l/L for critical point at A = 2. 
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FIG. 13. Effective conformal anomaly vs inverse system size l/L for various values of A 
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